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The dynamics of the Bose-Einstein condensate (BEC) in a double-well potential is of- 
ten investigated under the mean- field theory (MFT). This works successfully for large 
particle numbers with dynamical stability. But for dynamical instabilities, quantum cor- 
rections to the MFT becomes important [Phys.Rev.A 64, 013605(2001)]. Recently the 
adiabatic dynamics of the double-well BEC is investigated under the MFT in terms of 
a dark variable [Phys.Rev.A 81, 043621(2010)], which generalizes the adiabatic passage 
techniques in quantum optics to the nonlinear matter-wave case. We give a fully quan- 
tized version of it using second-quantization and introduce new correction terms from 
higher order interactions beyond the on-site interaction, which are interactions between 
the tunneling particle and the particle in the well and interactions between the tunneling 
particles. If only the on-site interaction is considered , this reduces to the usual two- mode 
BEC. 
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1. Introduction 

The dynamics of the Bose-Einstein condensate (BEC) is often investigated under the 
mean-field theory (MFT), which provides a classical field equation for the nonlinear 
matter wave and is the classical limit of the quantum theory in large number limit. 
MFT works successfully in predicting experimental results as the particle numbers 
of BEC is very large, so the classical limit captures the essence. In the vicinity of 
a mean-field dynamical instability, however, the quantum correction to the MFT 
becomes important ™ They provide accurate predictions for the dynamics by com- 
bining the mean-field with the fluctuations. As most investigations are done under 
the MFT, it is interesting to look at their behavior under quantum corrections. 

On the other hand, the analogy to quantum optics is noticed in investigating 
BEC, because of its macrosscopic coherence that allowes us to view it as a large 
atom. A recent example is the discovery of a dark variable for the double- well BEC 
by Ottaviani et al which is a nonlinear matter- wave extension of Vitanov and 
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Shore's work ^ of a two-state atom, where a dark variable is found by looking 
into the analogy between the optical two-state Bloch equation and the three-state 
stimulated Raman adiabatic passage (STRAP) equations 1 They compose three 
variables analogous to the Bloch spin from the ground states of the two wells and 
the dark variable is found by comparing the equation of the three variables with 
the STRAP equation. The adiabatic splitting, transport and self-trapping of the 
double- well BEC are investigated using this dark variable ^ under the MFT. The 
STRAP can be utilized to transport BEC in multiple wells EHHl 

In the paper, we give a fully quantized version of the work of Ottaviani et al. We 
introduce further quantum corrections by including higher order interactions beyond 
the on-site interaction. If these quantum corrections are turned off, th is mode l would 
reduce to the two-mode BEC widely used in quantum entanglement 11QH1SI 
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2. The on-site Interaction 

In MFT, the two-mode approximation is applied to the GP equation and the dy- 
namics of the double- well BEC ^ is 

where the Hamiltonian is given by 

U -\ n e R + U R \c R f)- 
Introducing the field operator ^ (r , t) , the Hamiltonian can be written as 

rl = Hq + rli n t 

+ |y d 3 r* + * t **, (3) 

where Ho = — V 2 /2m + V (f, t) is the single particle Hamiltonian with V (r, t) the 
double- well potential and g = 47ra s /m is the nonlinear interaction parameter with 
a s the s-wave scattering length. Only spatial degrees of freedom and two-body in- 
teractions are considered. Under the two mode approximation, where all modes 
are omitted except the condensate modes, the field operator can be expanded 
i& (r, t) = cll (t) 4>l (r) + a R (t) 4> R (r), where &l,r and 4>l,r are the annihilation op- 
erator and the ground state of the left and the right well respectively. The overlap 
between the two ground states is neglected because it is small compared with the 
on-site part, §* L 4> R \4>l\ ■ Now the second quantized Hamiltonian becomes 

TL = e L a' L a L + ila [ L a R + —a' L a' L a L a L 

+ e R a R a R + fla R a L + ^a R a R a R a R , (4) 
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with the parameters 

cl,r = J d 3 r (f>* L R H (j)L,R, (5) 

= J d 3 r 4>* L Ho(f>R = J d 3 r<t>* R H <l> L , (6) 

Uo = 9 I d 3 r |0 L | 4 = / d?r \^ R \\ (7) 



with cl,r the chemical potential of the two wells respectively and Uo the nonlinear 
on-site interaction between particles in the same well. The tunneling rate f2 be- 
tween the two wells is negtive, but it can also be defined to be positive by adding 
corresponding minus signs in the Hamiltonian. 

It is convenient to rewrite the Hamiltonian by the Schwinger angular mo- 
mentum operators J x — (a^ R ai, + o^clr) /2, J y = (a R az, — a^ORj /2i, J z = 

(^a R aii ~ o} L aj^j /2, where J x ^ y corresponds to the correlation between the two wells 

and J z is the particle number difference between the two wells. The Hamiltonian 
under these angular momentum operators becomes 

E « ~ / M N 2 N 

H = -N + eJ z + 2QJ X + U (j 2 z + — - — 

= UoJ 2 z + eJ z + 2QJ X , (8) 

with N = Nl + N R the total particle number, E = + e R the sum of the two 
chemical potentials and e = e R — the difference of the two chemical potential. 
Terms containing N and E are assumed to be conserved and are neglected from the 
Hamiltonian. This is the often used Hamiltonian in quantum entanglement, which 
only the on-site interaction is included. The dynamics of the system manifests in 
the evolution of these angular momentum operators 

d (J*\ ( -iU -(e + 2C/ J 2 ) \ (j x \ 

Jy = e + 2U J z -iU -2Q N» " ^ 
\Jz J \ 2Q, / \JzJ 

Equation "Eq. ([S])" is the second quantized version of the dynamics given in Ref. [2] 
where it is described by the evolution equation of the Bloch spin vector ^1 

-(e + Uu) 
(e + UuS) -2u) | ( v | , (10) 
2uj 



dt 



d 
dt 
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The elements of the Block spin vector fi, i/and uj are composed from the two modal 
populations. 

The two evolution equations would resemble each other with the angular momen- 
tum operator corresponding to the Bloch spin vector, if we normalize the angular 
momentum operator by N/2, 2J x>VtZ /N ~ u,v,w. The nonlinear interaction param- 
eter becomes g = ATrNa s /2m under this normalization, which becomes the same 
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g in "Eq. (1101) " . Without the normalization, the adiabatic splitting, transport and 
self-trapping can be given by the motion of the angular momentum operator on the 
Bloch sphere with radius j = N/2, while in Ref. 2 ; it corresponds to the motion 
of the Bloch spine in the unit sphere. The adiabatic transport corresponds to the 
variance of J z from —N/2 to then to N/2 and the self-trapping corresponds to its 
variance from —N/2 to then back to N/2. 

Except the resemblance of their appearance, these two equations are quite dif- 
ferent. The first one is the quantized version of the second one and it is the evolution 
equation for quantum operators rather than classical variables. We can approximate 
the operators with their expectation values, with the lowest-order approximation 
corresponding the MFT and the second-order approximation corresponding to the 
equation given by Ref. [TJ Also two diagonal terms of — iIIq emerge in the quantized 
version, which puts a phase factor to J x and J y . So the phases are also important 
in the dynamics. The ph ases and particle numbers can be investigated under the 
mean field approximation ^H, where they exhibit oscillations in the phase space. In 
the angular momentum space, we can choose the eigenstates of J z as the basis set 
{\jm)}. The state of the system at time t is given by |^ (t)) = exp {—iHt} \^o), 
where l^o) = \y~T~) 1S * ne initial state with the left well population. The adiabatic 
transport of the system then means the evolution from \^ ) to \^ (t)) = |yt)- 

3. Higher Order Interactions 

The above quantized version of Ref. [2] only includes the nonlinear effects up to the 
on-site interaction. It would be interesting to investigate the effects of the neglected 
interactions, which becomes useful when its behavior is investigated under a wide 
parameter regime. To do this, we add to the Hamiltonian "Eq. (U)" the overlapping 
part 4>* L (j}R an d introduce two new parameters Ut and Uu to describe their effects, 
where Ut captures interactions between the tunneling particle and the on-site par- 
ticle and U u captures interactions between the tunneling particles themselves. We 
term the tunneling particle as 'tunnelier' for simplicity. Now the full Hamiltonian 
is given by 

H = e L a [ L a L + ila' L a R + —a' L a' L a L a L 

(ft ft \ Uff /ff ft 

a' L a' L a L a R + a* L a R a L a,L J + — [a^ L aRaR + 2a , L a R anaL 

. -t - i o-t - ^0 -t -t - * 

+ £Ra R aR + Ma R a L —a R a R a R a R 

+ U t (a R a R a R a L + a R a) L a R a R ^j + (a R a R a L a L + 2a R a f L a L a R ^ , (11) 

with the newly introduced parameters given by Ut = gj d 3 r | 2 (4'l^r) ~ 
g ( d\ \^ R \ 2 and U tt = g ( d?r ^R? - g [ d 3 r (^l) 2 = 
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g J d 3 V tpltp Ripjtip l ■ For simplicity reasons, we have assumed ipj^ipR = i])\ij)L- In 
the Schwinger representation, this Hamiltonian becomes 

H=jN + eJ z + 2QJ X + U ^ + ^- - y 

+ 2Ut (N-l)J x + -2- (j-J- + J+J+ + 2 ( .7. J- + ./_ ./_ j - 2 A 

/ 7y2 at 

= —AT + eJ z + 2^4 + Z7 I J? + — - — 



iV 



2 • t,u \ v * , 4 2 

+ 2C/ t (N - 1) 4 + C/t t (j 2 - J 2 _ 2 J| + ^ 

= (U - 2U tt ) J 2 Z +eJ z +2(n + U t (N - 1)) J x + U tt (jj - J y 2 ) . (12) 

This representation illustrates the interactions between the tunneliers, with J_ J_ 
the interaction between left tunneliers (particles tunneling from the right well to 
the left well), J+J+ the interaction between the right tunneliers (particles tunneling 
from the left well to the right well) , and J+ J_ + J_ J+ the interaction between the 
left tunnelier and the right tunnelier. We can see the modifications introduced to the 
Hamiltonian by comparing it with "Eq. (|SJ)" . The interaction between the tunnelier 
and the on-site particle adds a term Ut (N — 1) to the original tunneling parameter 
fi, the interaction between the left and the right tunneliers adds a term —2Uu 
to the original nonlinear on-site interaction parameter Uq and the interaction of 
the left tunneliers and that of the right tunneliers add another nonlinear terms 
Uu {jx — Jy^j to the whole Hamiltonian. The equation of motion of these operators 



is 



d 
dt 



( i x 

Jy 



1 -i(U -Utt) -(e + 2(U a -Uu)J z ) \ lj x 

e + 2 (U - 3U tt ) J z -i (U ~ 3U tt ) -2 (Q + U t (N- 1)) J y 

y 2(n + U t (N-l)+2UttJx) -iZU« ) \J* 

(13) 

Compared with "Eq. @" , this equation is non-symmetric and there is even a 
phase term of J z . By the definition of J 2 , there is no phase term as the phases 
carried by d^ is canceled out in a^et£ and d R aR. The origin of this phase term of 
J z manifests in the commutation relationship of J x and J y , which is introduced in 
deriving this equation of motion from the Hamiltonian "Eq. (fT2"j)" . 

In Ref. [21 the nonlinear parameter Uq generates a rich adiabatic dynamics of the 
double- well BEC. Here the quantum corrections introduced by Uu and Ut makes the 
dynamics even more complex. When Uu and Ut are small compared with Uo, the 
influence may not be very explicit; but when they are large enough, the formulations 
above would lose sense as the two-mode assumption no longer holds true. Numerical 
simulations are needed to investigate their influence to the dynamics under various 
parameter regimes. 
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4. Conclusion 

We have obtained quantum corrections of higher order interactions for the quan- 
tized dynamics of the double-well BEC and given a fully quantized expression of it 
under the two-mode approximation. Two new parameters are obtained to express 
the influence of the interactions between the tunneling particle and the on-site par- 
ticle and between the tunneling particles. So the dynamics of the double-well BEC 
becomes more rich under this quantized version. This allowes numerical simulation 
for a wide range of parameter regime and new phenomenon may be obtained. 
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